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SUMMARY: It has recently been reported that some of the flux density val- 
ues for an evolved supernova remnant (SNR) HB 3 (G132.7+1.3) are not accurate. 
In this work we revised an analysis of the possible thermal emission at radio fre- 
quencies from the SNR HB 3 using the recently published, corrected, flux density 
values. A model including a sum of non-thermal (purely synchrotron) and ther- 
mal (bremsstrahlung) component is applied for fltting integrated radio spectrum 
of the SNR. The contribution of thermal component in total volume emissivity at 
1 GHz was estimated to be ~ 37%. The ambient density was also estimated to 
be n « 9 cm"3 for the T = 10"' K. Again, we obtained the relatively signiflcant 
presence of thermal emission at radio frequencies from the SNR so we could sup- 
port interaction between SNR HB 3 and adjacent molecular cloud associated with 
the H II region W3. Our model estimates for thermal component contribution to 
total volume emissivity at 1 GHz and ambient density are similar to those obtained 
earlier (~ 40 %, ~ 10 Cm It is clear that the corrected flux density values do 

not change the basic conclusions. 
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1. INTRODUCTION 



The presence of thermal emission at radio fre- 
quencies may be a useful tool for identifying in- 
teractions between supernova remnants (SNRs) and 
molecular clouds, and also for estimating the ambi- 
ent density near SNRs using radio continuum data 
(Urosevic & Pannuti 2005, Urosevic, Pannuti & 
Leahy 2007). In this work we argue for the pres- 
ence of a thermal bremsstrahlung component in the 
radio emission from SNR HB 3 in addition to the 
synchrotron component. It is possible that SNRs can 



be sources with significant amount of thermal radia- 
tion and as Urosevic & Pannuti (2005) stated, there 
are two basic criteria for the production of a sig- 
nificant amount of radio emission through the ther- 
mal bremsstrahlung process from an SNR: the SNR 
evolves in an environment denser than the average 
and its temperature must be lower than the average 
(but always greater then the recombination temper- 
ature). Two cases have been considered: thermal 
emission at radio frequencies from a relatively young 
SNR evolving in a dense molecular cloud environ- 
ment (n w 100 — 1000 cm~^) and extremely evolved 
SNR (approximately 10^ — 10^ years old) expanding 
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in a dense warm medium (n « 1 — 10 cm '^). For 
a detailed discussion about the issue see Urosevic 
& Pannuti (2005) and Urosevic, Pannuti & Leahy 
(2007). 

Urosevic, Pannuti & Leahy (2007) have ana- 
lyzed the broadband (22 - 3900 MHz) radio spec- 
trum of the Galactic SNR HB 3 and have discussed 
the possible thermal radio emission from the SNR. 
Earlier published observations have revealed that a 
curvature is present in the radio spectrum of SNR 
HB 3, indicating that a single synchrotron compo- 
nent appears insufficient to adequately fit the radio 
spectrum. Urosevic, Pannuti & Leahy (2007) have 
suggested that more natural explanation for the ap- 
parent spectral index variations found by Tian fc 
Leahy (2005) is synchrotron emission, which domi- 
nates at lower frequencies, and bremsstrahlung emis- 
sion, which dominates at higher frequencies. They 
have found w 40 % for the thermal component con- 
tribution to total volume emissivity at 1 GHz and 
have also estimated the ambient density, implied by 
the presence of thermal component, to be ~ 10 cm~^. 
Recently, Green (2007) has reviewed the radio spec- 
trum of SNR HB 3, noting the difficulties in deriving 
accurate flux densities for the remnant, particularly 
at higher frequencies, due to thermal emission from 
the nearby bright Hii region W3 (IC 1795) and its 
surroundings. He pointed out that some of the ear- 
lier published flux density values used by Urosevic, 
Parmuti & Leahy (2007) in their analysis are not ac- 
curate and that the spectrum of the SNR is well rep- 
resented by a simple power low spectrum as well as 
the contamination with thermal emission from ad- 
jacent regions is the cause for the reported spectral 
flattening of the spectrum. In this work we present 
the results of our analysis using the corrected, re- 
cently reported by Green (2007), flux densities for 
SNR HB 3. 



2. The model 



A spectrum of SNRs in radio domain is usu- 
ally represented by an ordinary power law. If the 
frequency is in GHz, the flux density can be repre- 
sented by the following expression: 

Sv = Sigut. ■ (1) 

where: S'ighz is the ffiix density at 1 GHz, and a 
is the radio spectral index. In order to distinguish 



the contribution of thermal and non-thermal com- 
ponent within the total radiation, SNR radio inte- 
grated spectrum was fitted by a simple sum of these 
two components. If the frequencies are still in GHz, 
the relation for the flux density can be written as 
follows: 

c-T 

o _ oNT / -a , '^IGHz ,.-0.1^ fr,^ 
— '-'iGHz ['^ + cNT ^ i' 
'^IGHz 

where: SJq^^ and Sj^^jj^ are flux densities corre- 
sponding to thermal and non-thermal component, 
respectively. The spectral index is considered to be 
constant in the SNR shell. It is also considered that 
the thermal radiation is optically thin and has the 
spectral index equal 0.1 at any point. As the radio 
frequency increases, the amount of synchrotron ra- 
diation from an SNR decreases and the amount of 
thermal bremsstrahlung emission becomes more sig- 
nificant. In our model it is also considered that the 
synchrotron radiation, optically thin at any point, is 
not absorbed or scattered by thermal gas. 

The volume emissivity of thermal 
bremsstrahlung radiation for an ionized gas cloud 
is proportional to the square of the electron (or ion) 
volume density n, i.e.: 

= 7 X 10-^VT-5, (3) 

where n is in cm~^ and thermodynamic temperature 
T in K. Having determined total e,^ and thermal com- 
ponent contribution to total volume emissivity, the 
density of the of the interstellar medium (ISM) can 
be estimated using equation 3. 

This model is valid only in the approximation 
of constant density and temperature. The model it- 
self also presumes a simple sum of non-thermal and 
thermal component, while the attention is not paid 
to the fact that the dependence of flux density could 
be another, more complicated, function of thermal 
and non-thermal components. It is also important 
to note that this model, in general, does not dis- 
tinguish between thermal and non-thermal emission 
with the same spectral index (i.e. the case of lower 
synchrotron spectral index). 

Despite these drawbacks, our model represents 
a useful tool for estimating the contribution of ther- 
mal bremsstrahlung component to the total volume 
emissivity and ambient density using radio contin- 
uum data. 
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3. SNR HB 3 (G132.7+1.3) 



From the Green (2006) paper, the value of the 
radio spectral index is 0.4, Siguz = 45 Jy, size is 
around 80 arcmin and the SNR is S (shell) type. 
On the other hand, a combination of radio shell 
morphology with a center-filled thermal X-ray mor- 
phology has led to the classification of SNR HB 3 
as a mixed- morphology SNR (Urosevic, Pannuti & 
Leahy 2007 and references therein). It is the one 
of the largest SNRs currently known (Fcscn et al. 
1995, Tian & Leahy 2005). The distance from the 
SNR is about 2.2 kpc (Tian & Leahy 2005, Shi 
et al. 2008). SNR HB 3 size, based on a dis- 
tance of 2 kpc is 60 X 80 pc (Tian & Leahy 2005 
and references therein). Lazendic & Slane (2006) 
stated that the SNR is 90 x 120 arcmin in diam- 
eter. Kovalenko, Pynzar & Udal'tsov (1994) have 

reported: a = (0.51 ± 0.12), Fcscn ct al. (1995): 
a = (0.64 ± 0.01) (also pointed by Lazendic & Slane 
2006) and Landecker et al. (1987): a = (0.60±0.04). 
Green (2007) has found a = (0.56 ± 0.03). Tian & 
Leahy (2005) pointed spectral index variations with 
most values between 0.3 and 0.7. 

Shi et al. (2008) have extracted 4800 MHz 
total intensity and polarization data of HB 3 from 
the Sino-German 6 cm polarization survey of the 
Galactic plane made with the Urumqi 25 m tele- 
scope, but they could not give a total flux density 
at 4800 MHz of the whole SNR because of a low res- 
olution. They have found a radio spectral index of 
HB 3 of a = —(0.61 ± 0.06) using only three flux 
densities, at 1408 MHz, 2695 MHz and 4800 MHz, 
and have concluded that there is no spectral flatten- 
ing at high frequencies. Shi et al. (2008) have also 
pointed out that a reliable observations of SNR HB 3 
at frequencies above 3000 MHz are crucial to confirm 
a spectral flattening. 

A radio pulsar PSR J0215-h6218 has been dis- 
covered within (in projection) the SNR HB 3 bound- 
ary, but it appears to be much older than the rem- 
nant and therefore not associated with the SNR 
(Lazendic & Slane 2006, Lorimer, Lyne & Camilo 
1998). 

Fesen et al. (1995) pointed that SNR HB 3 
is relatively opticaUy faint SNR. Difi^use and fila- 
mentary optical emission has been detected from the 

SNR, with the strongest emission along the western 
SNR sheU (Lazendic & Slane 2006). Reich, Zhang 
& Fiirst (2003) pointed that the SNR radio sheU is 
located at the western edge of the H ii region com- 
plex W3-W4-W5. Optical emission from the SNR 
was found to be well-correlated with the radio emis- 
sion, with a multiple shock structure found in the 
western SNR shell and lack of emission in the south- 
east region (Lazendic & Slane 2006 and references 
therein) . 

Most of the mixed-morphology SNRs are in- 



teracting with molecular or H I clouds, as indicated 
in some cases by infrared line emission or OH masers 
(Rho & Petre 1998). OH (1720 MHz) masers, which 
are recognized diagnostic for a molecular cloud 
interaction with a SNR, have been detected towards 
the W3/HB 3 complex (Lazendic & Slane 2006, Ko- 
ralesky et al. 1998). 

Urosevic, Pannuti & Leahy (2007) pointed 
that the X-ray emission is seen to lie entirely within 
the radio shell of HB 3. To obtain an independent es- 
timate of the ambient density of the ISM surrounding 
HB 3, they have performed spectral fitting on the ex- 
tracted j4SCA GIS spectra and have calculated elec- 
tron densities to be rig « OA f"^^^ cm'^ for the cen- 
tral region and rie « 0.1/~^/^ cm~'^ for the northern 
and southern regions (see Table 2 in Urosevic, Pan- 
nuti & Leahy 2007), where / represents the volume 
filling factor. 



Analysis and results 



Green (2007) pointed out that the radio ob- 
servations of the SNR HB 3 are complicated by con- 
fusing thermal emission from an adjacent Hii region 
W3. He has shown that some of the earlier pub- 
lished fiux density values (used by Urosevic, Pannuti 
& Leahy 2007 in their analysis) are not accurate. 
Green (2007) also pointed the huge problem of deriv- 
ing the accurate flux density values and has listed the 
corrected values for SNR HB 3. Particulary, he re- 
ported corrected values for 408 and 1420 MHz points 
and also corrected uncertainty for 865 MHz point, all 
from Tian & Leahy (2005) and references therein. He 
also derived an integrated flux density for the Effels- 
berg 2695 MHz survey data. Green (2007) excluded 
3650 and 3900 MHz points (from Tian & Leahy 2005 
and references therein) due to their possible contam- 
ination with thermal emission associated with W3. 
He also rescalcd 22, 38 and 178 MHz points, used 
by Urosevic, Pannuti & Leahy (2007) in their anal- 
ysis, to be on the scale of Baars et al. (1977). We 
have revised an analysis of the possible thermal emis- 
sion contribution in total volume emissivity at radio 
frequencies from an evolved SNR HB 3 using data 
points for integrated radio flux density from Green 
(2007) for a range of 22 MHz to 2.695 GHz (see Table 
1 in Green 2007). 

The parameters of our model fit (equation 2) 
are as follows: a = (0.70 ± 0.05), S^^^^ = (30.45 ± 

5.22) Jy, 5?GH./5rGHz = (0.59±0.27),x?ed = 0-19, 9 
degrees of freedom (dof) as it can be seen in Table 1. 
Xred represents reduced (x^/dof). The parame- 
ters of purely non-thermal model fit (equation 1) are 
as follows: a = (0.56 ±0.02), Sf^^^^ = (49.18 ±2.06) 
Jy, Xred ~ 0.31, 10 dof as it can be seen in Table 2. 
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Table 1. The fit paramc't(TS for our model 
for SNR IIP.;-!. 

0.70 ±0.05 30.45 ±5.22 0.59 ±"6.27 0.19 (9)~ 



Table 2. The fit parameters for purely non- 
thermal model for SNR HB3. 

a ^^f^Hz (Jy) xLi (dof) 
0.56 ±0.02 49.18 ±2.06 0.31 (10) 



It can be noted that the radio spectral index 
value is higher than the value from Green (2006) both 
for the purely non-thermal {a = 0.56 ±0.02) and our 
model fit (a = 0.70 ± 0.05) calculations. The results 
from the purely non-thermal model fit responds to 
the values placed in Green (2007). Our model fit 
radio spectral index estimate is closer to the value 
placed in Lazendic & Slane (2006) and Fesen ct al. 
(1995) than in Green (2006, 2007). Figure 1 shows in 
a solid line fit by non-thermal plus thermal model, 
while the dotted line shows the fit by purely non- 
thermal model through the given sample. 



the relatively significant presence of thermal emis- 
sion at radio frequencies from the SNR. Our model 
estimates for thermal component contribution to to- 
tal volume emissivity at 1 GHz (« 37 %) and am- 
bient density (« 9 cm~'^) are similar to those ob- 
tained earlier (rs 40 %, « 10 cm~^) by Urosevic, 
Pannuti & Leahy (2007). The fact that essentially 
the same thermal component again minimizes the ) 
indicates that its presence cannot be ruled out by the 
corrections to the flux densities. 

If we assume the value for the compression 
parameter to be 4 we can roughly estimate pre- 
shock ISM number density as no « 2.25 cm^"^ for 
T = 10" K. 

It is clearly visible that our ambient density es- 
timates could support the possibility that the SNR is 
indeed expanding in a dense ISM. Based on our anal- 
ysis we can support the possibility that SNR HB 3 is 
indeed interacting with the molecular cloud material. 
The presence of the thermal bremsstrahlung compo- 
nent in the radio spectrum of SNR HB 3 suggests 
that this SNR is in fact interacting with adjacent 
molecular cloud associated with the Hii region W3. 

It should be pointed out that the further mea- 
surements at the highest radio frequencies (> 3 GHz) 
are required for a detailed discussion of the issue. 




v(GHz) 



Fig. 1. The integrated spectrum of HB 3. The 

full line represents the fit by non-thermal plus ther- 
mal model, while the dotted line represents fitting by 
purely non-thermal model. 



If a mean value of 37 %, from our analy- 
sis, for the thermal component contribution to to- 
tal volume emissivity at 1 GHz is assumed, we get 

n « 9 cm^'' for the assumed post-shock tempera- 
ture value of 10^ K. We have assumed d = 2 kpc 
and D = 70 pc for consistency with Urosevic, Pan- 
nuti & Leahy (2007). It is clearly seen that for the 
corrected integrated flux density values we also get 



Conclusions 



In this work we revised an analysis of the pos- 
sible thermal emission at radio frequencies from an 
evolved SNR HB 3. Some of the earlier published flux 
density values for SNR HB 3 are shown to be non- 
accurate. Here we present the results of our analysis 
using the recently published, corrected, flux densities 
so the main conclusions are: 

1. The contribution of thermal component in to- 
tal volume emissivity was estimated to be 
w 37% and the ambient density was also esti- 
mated to be n w 9 cm^"^ for the T = 10** K. 

2. Our model estimates for thermal component 
contribution to total volume emissivity at 
1 GHz and ambient density are similar to 
those obtained earlier (« 40 %, « 10 cm~^). 
It is clear that the corrected flux density val- 
ues do not change the basic conclusions. 

3. The presence of the thermal bremsstrahlung 
component in the radio spectrum of SNR HB 
3 suggests that this SNR is in fact interacting 
with adjacent molecular cloud associated with 
the Hll region W3. The presence of thermal 
emission at radio frequencies may be a useful 
tool for identifying interactions between SNRs 
and molecular clouds and also for estimating 
the ambient density near SNRs using the radio 
continuum data. 

4. The lack of data at higher radio frequencies 
unable us to give a flrmer conclusion about 
the issue. 
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AHAJIH3A MOryiiE TEPMAJIHE EMHCHJE HA 
PAjmO *PEKBEHHMJAMA EBOJIYHPAHOr OCTATKA XB3 

Jl. OHHti^ H XL ypomeBnti^ 

^ Kamedpa aa acmpoHOMujy, MameMamuuKu (fiaKyAmem, yHueepaumem y Eeozpady 
CmydeHmcKU mpz 16, 11000 Eeozpad, Cp6uja 

UDK 

IIpemxodHO caonmmeibe 



He/iaBHO je o5jaBJi.eHO ^a cy neKe Bpe/i;- 
HOCTH rycTHHa (jjJiyKca eBOJiynpaHor ocTaiKa cy- 
nepHOBC (OCH) HB 3 (G132.7+1.3) HCTa^HC. V 
OBOM pa^y CMO hohobhjih anajiHsy Morytie lep- 
Majine pa/iHO-eMHCHje OCH HB 3 KopuciehH 
He^aBHO ny5jiHKOBaHe, KopHroBane Bpe^HOCTH, 
rycTHHa (|)JiyKca. Mo^eji KojH npeinocTaBJta 
cyMy HeiepMajiHe h TepMajine KOMnoneHTe je 
npHMeaeH 3a ^JHTOBaite pa/iHO-cneKTpa ocTaiKa. 
IIpHcycTBO TepMaJiHB KOMnoHCHTC y yKynHoj sa- 
npcMHHCKoj eMHCHBHOCTH Ha 1 GHz je npoAeHjeno 



Ha « 37 %. FycTHHa okojihc cpejxyme je TaKoije 
npon;cHjHa na n « 9 cm^'^ sa T = 10* K. Hohobo 
je o^pcI)CHO SHaMajHO npncycTBO TcpMajiHC kom- 
noHCHTC y yKynHoj aanpcMHHCKoj cmhchbhocth 
TaKO Ji& MO»ccMO jja noziynpcMO XHnoTcay o 
HHTepaKn,iijii HSMcljy ocTaTKa h MOJicKyjiapnor 
o6jiaKa. IIpon,CHC npHcycTBa TcpMajiHc komho- 
HBHTe y yKynnoj sanpeMHHCKoj cmhchbhocth na 
1 GHz H rycTHHC okojihc cpc^HHC cy cjih^hb ca 
paHHje o/npe^eHHM (« 40 %, « 10 cm~^). Jacno 
je BH^jbHBO fla Bpe^HOCTH KopHroBaHHx rycTHna 
(|)jiyKca He Meaajy ocHOBne saKJBy^iKe. 
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